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Abstract 

We study the energy level statistics of the states in E(5) and X(5) dynamical 
symmetries. The calculated results indicate that the statistics of E(5) symmetry 
is regular and follows Poisson statistics, while that of X(5) symmetry involves two 
maxima in the nearest neighbor level spacing distribution P{s) and the A 3 statistics 
follows the GOE statistics. It provides an evidence that the X(5) symmetry is at 
the critical point exhibiting competing degrees of freedom. 
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The nature of “shape” phase transitions in finite quantal system is a fundamental 
issue and has been the subject of many investigations. In the last years, there have been 
many preeminent works concerning this subject in nucleus, including both theoretical 
calculations |l|, H, ^ and experimental results^, ||. Particularly, two symmetries at 
the critical points of shape phase transition, called E(5) and X(5), have recently been 
proposed]^, It is a major breakthrough in the study of critical point behavior of 
nucleus undergoing a shape phase transition because it gives us simple analytic solutions 
while historically we have to resort to numerical calculations to describe such nucleus. 

Empirical examples of such kind of symmetries have been found soon after they 
were put forward. For instances, the ^°®Pd, Ru isotopes are proposed 

to be the ones in E(5) symmetry[y, [^, ^ and the ^®^Sm, ^®°Nd, ^°^Mo in X(5) 

symmetry |]I^, |^. These examples show that the obtained new symmetries represent 

a helpful theoretical tool which could well describe the structure of the realistic nucleus 
at the critical point of a shape transition, in a sense complementing the role played by 
the dynamical symmetries of the Interacting Boson Model (IBM) [l^, However, there 

remain some interesting aspects on the statistical properties of the E(5) and X(5) sym¬ 
metries. (1) Whether the nuclear system displays ordered or chaotic spectrum in the new 
symmetries? (2) Does the statistics give a strong evidence to support that the obtained 
new symmetries can well describe the nucleus at the critical point in a shape transition? 
(3) The new symmetries are based on the particular potential amenable to analytic de¬ 
scriptions which well approximate to the “true” potential. Does the simplification of the 
potential eliminates part of the statistical properties of the realistic nucleus at the critical 
point? Or more concretely, are the results consistent with the statistical properties ob¬ 
tained by Alhassid and collaborates)^, in the framework of IBM? Therefore it is 

crucially important to obtain the statistics of the new symmetries and answer the above 
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opening questions. In this sense, we explore the statistical properties of E(5) and X(5) 
symmetries in the paper. 


H 


We start from the Bohr Hamiltonian [pT 
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In the E(5) symmetry, the potential V{f3, 7 ) depends only on the deformation parameter 
[3 and it can be written as I7(/9, 7 ) = f/(/3). Considering the case, in which the potential 
is simplified as an infinite well, one has the eigenvalues P 
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where H is a parameter, Xi,T is the .^th zero of the Bessel function Jp{z), and the order 
of the Bessel function is z/ = r + | with r being the quantum number similar to the 
irreducible representation of 0(5) group in the IBM, /3w is the width of the well. 

In the X(5) symmetry, the potential was supposed to be a square well in the variable 
(3 and a harmonic oscillator in 7 with no /3 — 7 couplings, which is a good approximation 
in the U(5)-SU(3) shape transition. In this case, the eigenvalues are given as[§] 


E{s, L, K, M) = ^ 


(3) 


where 



(4) 


and 



{n^ + 1 ) 


4 {K/2f 
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The Xs,L is the sth zero of Bessel function Ju{z) of (irrational) order z/ = [J(J + l)/3 + 
9/4]02, L is the total angular momentum (with projections K on the symmetry axis and 
M on the quantization axis in the laboratory frame), a is a parameter related to the 
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interaction strength of the potential = ( 3 a)^ 7^/2 and is the number of 7 -vibration 
quanta, while /9w is the same as that in E(5) symmetry and the (/3^) is the average 
of over the wave functions. Combining all variables, one obtains the most general 
expression [|] 

E(s, L, n^,K,M) = Eo + B{Xs,l)‘^ + + CK"^ , ( 6 ) 

where 



Before analyzing the energy level statistics, we estimate the relative value of the param¬ 
eters A/B m the light of the “true” potential given in the Fig.l. of Ref.[§. We can see 
from the hgure that the width of the inhnite well to approximate the “true” potential is 
/3w ~ 0.8, and most of the quantum states can thus be restricted in Uj^ < 0.5 (arbitrary 
unit). If we assume that the total energy of the quantum states are conserved m. E = 0.5, 
the quantum states of the 7 degree are then mainly restricted in < 0.5, as the potential 
is completely decoupled into two separate parts in f3 and 7 respectively. Since the X(5) 
symmetry is obtained around 7 = 0 |p], the corresponding restriction in 7 by the total 
energy must be very small. Taking this small value of 7 as 0.1, we then work out that the 
strength of the two dimensional oscillator (3a) must be more than 10 according to the 
formulation = ( 3 a)^ 7 ^/ 2 . With Eqs.(7) and ( 8 ), we obtain that the parameter A must 
be more than 25 if we take as 0.4 and the parameter B = 1.56. Consequently we 

get the relative value A/B > 16. 

Our estimate is consistent with the recent empirical example of X(5) sym¬ 

metry, where the nucleus presents the X(5) pattern not only in its ground state band, but 
also in its low-lying = 1,2] K = 2n^ bands. From the level scheme given in the Fig.2. 
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of Ref.||15||, we can obtain the empirical values of the parameters, A ^ 800 and B ^ 30 
(in keV). The relative value is clearly in the region of A/R > 16 as we have proposed. 
Such a consistence supports that the assumption to obtain X(5) symmetry around 7 = 0 
in the Ref.|^ is quite reasonable since in fact, most of the quantum states in the realistic 
nucleus at the critical point of the transition from spherical to axially-deformed shape are 
restricted in 7 < 0 . 1 . 

To analyze the energy level statistics in the dynamical symmetries, we take the fol¬ 
lowing process. At hrst with Eqs.(2) and ( 6 ), we calculate the energy levels of the nucleus 
in E(5) and X(5) symmetries, respectively. The parameter ^ and r in E(5) symmetry 
are truncated manually at a maximal number and We also truncated s and 

at and respectively in X(5) symmetry. After the spectrum {Ei} has been 
determined, it is necessary to separate the smoothed average part whose behavior is 
nonuniversal and can not be described by random-matrix theory(RMT) . To do so we 
count the number of the levels below E so that one can dehne a staircase function N{E) 


of the spectrum (see for example Ref. |]lq ) 


N{E) = Na,{E) + . (10) 

Then we £x the Nav{Ei) semiclassically by taking a smooth polynomial function to £t the 
staircase function N{E). We obtain finally the unfolded spectrum with the mapping 

{E,} = N{E,). (11) 

This unfolded level sequence {E^} is obviously dimensionless and has a constant average 
spacing of 1 , but the actual spacings exhibit frequently strong fluctuation. 

We have used two statistical measures to determine the fluctuation properties of 
the unfolded levels: the nearest neighbor level spacings distribution(NSD) P{S) and the 
spectral rigidity A^^L). The nearest neighbor level spacing is dehned as Si = (Ej+i 
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The distribution P{S) is defined as that P{S)dS is the probability for the Si to lie within 
the infinitesimal interval [S, S + dS\. For a regular system, it is expected to behave like 
the Poisson statistics 




( 12 ) 


whereas if the system is chaotic, one expects to obtain the Wigner distribution 


P(S) = (ir/2)S exp(-irS'V4), 


(13) 


which is consistent with the GOE statistics [p^, 
the Brody distribution 


With the Brody parameter u in 


P^(S) = a(l + o;)^“exp(-a^(^+“)), 


(14) 


where 


a = T[{2 + uj)/{l + uj)Y/^ (15) 

and r[x] is the P function, the transition from regularity to chaos can be measured with 
the Brody parameter oj quantitatively. It is evident that a; = 1 corresponds to the GOE 
distribution, while a; = 0 to the Poisson distribution. A value 0 < a; < 1 means an 
interplay between the regular and the chaotic. 

As to the spectral rigidity A^^L), it is defined as 

A3{L) = (^minA,BY J ~ ~ ■, (16) 

where N{x) is the staircase function of a unfolded spectrum in the interval [—L/2,x]. 
The minimum is taken with respect to the parameters A and B. The average denoted 
by (■ ■ ■) is taken over a suitable energy interval over x. Thus from this definition A^lL) 
is the local average least square deviation of the staircase function N{x) from the best 
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fitting straight line. For the GOE the expected valne of A^i^L) can only be evalnated 
nnmerically, bnt it approaches the valne 

As{L) ^ {InL - 0.0687)/7r2 (17) 

for large L. And for Poisson statistics 


A3(L) = L/15. 


( 18 ) 


In onr calcnlation of E(5) symmetry, all the degenerate states cansed by additional 
qnantnm nnmber are considered jnst as one single state[^. That is becanse we mann- 
ally introdnced the additional qnantnm nnmber to distingnish the degenerate states with 
the same quantum number r. The degenerate states with different additional quantum 
numbers are statistically uncorrelated since no interaction is involved to link such states. 
If we do not remove the degeneracy, we would take the mixed ensembles (with different 
additional quantum numbers which are good quantum numbers in the dynamical symme¬ 
tries) into consideration and obtain the over-Poisson distribution in practical calculations 
since nearly 1/4 of the spacings of all are zero. In order to evaluate the fluctuations of the 
energy levels in one pure ensemble, these large amount of the statistically uncorrelated 
states should be removed. 

At hrst, we analyze the energy levels given in Eqs.(2) for the E(5) dynamical symme¬ 
try. The numerical results of the NSD and A 3 statistics show that the statistical feature 
does not depend on the angular momentum obviously. We then illustrate only those of 
the states with low spin-parity J'^ = O’*" but different sets of truncated numbers and 
t-M of E(5) symmetry in Fig.l (a), (b) and (c). Meanwhile the Brody parameter ui of 
the NSDp^ is also evaluated and listed in the hgures. The numerical results show that 
the statistics of the E(5) symmetry is quite close to Poisson statistics. When the number 
of the energy levels involved in the statistics (denoted by and increases, the 
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statistics becomes almost exactly the Poisson-type. It indicates that the E(5) symmetry 
is a dynamical symmetry which is completely integrable in classical limit. 


We then analyze the statistics in the X(5) symmetry. We can easily find that the 
energy eigenvalnes in the X(5) symmetry contain two terms, one is the sqnare of the zeros 
of the Bessel fnnction in f3, another is the solntion of two dimensional oscillator fnnction 
in 7. In order to display how the relative strength of the two terms affects the energy 
level statistics, we calcnlate the statistics in different valnes of the relative strength A/B, 
the resnlts are displayed in Fig.2. All of the relative strength A/B satisfy the condition 
A/B > 16 as we estimated above. It is apparent that when the parameter A/B is not very 
large, two maxima [p5[| appear in the NSD statistics (see the left part of Fig.2. (a), (b)), 
which is rather different from the finctnation properties in other cases except for those in 
the U(5) dynamical symmetry with collective backbending [p5|, [26[| . In theoretical point 
of view, the X(5) symmetry describes the critical point nnclei in a spherical to axially 
deformed shape phase transition analytically. The structnre of its energy levels must be 
not the same everywhere. In fact, the term B{xs,l)‘^ is approximately a rotational term 
of s, whereas the term An^ is a vibrational term of n^. The simultaneous appearance of 
these two terms induces a competition and makes the ensemble with the same angular 
momentum L and its projection on the symmetry axis K involves two sequences, one of 
which is in vibration, another is in rotation. When we unfold the above spectrum {Bj(L)}, 
the two sequences are normalized with an unique total average spacing, and their maxima 
in the spacing distributions are no longer the same, consequently two maxima emerge 
because of the uniform unfolding procedure. The same mechanism and results can be 
found in U(5) dynamical symmetry with collective backbending, where the the structure 
of the yrast band changes from the U(5) vibrational states to the rotational states with 
full d-boson number = N, which implies that a transition of collective motion mode 






may happen as the angular momentum increases. In this statistical point of view, the 
X(5) symmetry describes well the critical point behavior of nucleus undergoing a shape 
transition from U(5) dynamical symmetry to SU(3) dynamical symmetry which involves 
competing degrees of freedom. 

In the A 3 statistics when A is not very large (see the right part of Fig.2. (a), (b)), it 
is clearly that the statistics follows the GOE statistics, which suggests that the competing 
strengthes in the (3 and 7 degrees can cause the onset of chaos in X(5) symmetry. The 
case that the degree of chaoticity reaches its maximal when different competing strengthes 
coexist and are comparable can be also found in other works (see for example Ref. p5|, ^). 

Looking over Fig.2.(c) and (d), one can recognize that when we increase the strength 
of the oscillator in 7 degree of freedom, which corresponds to the parameter A, the 
two maxima in the NSD statistics soon vanishes and both the NSD and A 3 statistics 
become close to the Poisson statistics. It indicates that only the vibrational mode plays 
a dominant role in the system and the competition mentioned above does not exist. If we 
go to extremes to the case when the vibrational strength 3a —>■ cx), the 7 degree of freedom 
no longer exists. The X(5) symmetry becomes the E(5) symmetry in a four-dimensional 
space, and its statistics is certainly regular. 


The results obtained above are consistent with Alhassid and collaborates’ work 19 


The E(5) symmetry lies in the transitional region between a vibration to a 7 -unstable 
rotation, which is regular both classically and quantum mechanically in the work. It has 
been known that the Hamiltonian in the consistent Q formahsm|P 8 || of the IBM can be 
given as 

H = erid - kQ^(x) ' Q^ix), (9) 


where 


Qlix) = + Ss)l + x{d^d) 
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By changing e /k under restriction x = —/2, one can describe the U(5)-SU(3) transition 
with this Hamiltonian. To discuss the shape phase transition and relate it with X(5) 
symmetry, we reparametrize the ration e/n as (1 — C)/C- From the Ref. |l| we know that 
the flat potential near the critical point could be obtained when the parameter ^ G (0.025, 
0.029) [m in case that the boson number N = 10. We then transfer the parameter ^ into 
parameter t] used in Alhassid and collaborates’ work[^. The relation is 

1-C 

1-C + iVC’ 

and the value of rj near the critical point is 0.770 ~ 0.796. It is clearly that the critical 
point which corresponds to the X(5) symmetry lies in the most chaotic region(see the 
Fig.l. in Refill^). 

Finally, we turn to the question we have raised in the beginning. In X(5) symmetry the 
nucleus exhibits chaotic spectrum, which is probably caused by the comparable competing 
strengthes in the (3 and 7 degrees. The two maxima appearing in the NSD statistics 
suggest that the spectrum splits into two different sequences, one is rotational and the 
other is vibrational. Such a behavior illustrates that the nucleus in X(5) symmetry is at the 
critical point of the transition from spherical to axially-deformed shape with competing 
degrees of freedom. While the statistics in E(5) symmetry provides no evidence for such 
a drastic competition in the critical shape. That might be because the E(5) symmetry is 
a critical symmetry corresponding to an isolated point of second order phase transition. 
Very recently, it has been pointed out that such a second-order phase transition located 
at the triple-point of the three phases [^. The shape coexistence region shrinks as the 
parameter y in the Casten’s triangle changes from —\/l/2 to 0. In deed, as a general rule, 
phase coexistence can be found only in hrst-order phase transitions^]. There is no shape 
coexistence in E(5) symmetry and the spectrum does not split into two sequences, then 
no competition exists and its statistics is uniform(no two maxima in the NSD statistics). 
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The common subalgebra 0(5) in the U(5)—0(6) transition makes the statistics regular 
in this region and the ordered spectrum of nucleus in E(5) symmetry conhrms the past 
workJT^. It is worth to mention that the bifurcation effect (number of maxima more than 
one) in the NSD statistics might be a useful tool to determine whether the system is 
at the critical point of a hrst order shape transition or shape coexistence which involves 
competing degrees of freedom. 

In summary, we have calculated the statistics of E(5) and X(5) symmetries. The 
statistics of X(5) symmetry involves two maxima in the NSD statistics and the A 3 statis¬ 
tics follows the GOE statistics, which provides an evidence that the X(5) symmetry is at 
the critical point exhibits competing degrees of freedom. The statistics of E(5) symmetry 
is regular. That is probably caused by the common subgroup 0(5) in the U(5)—0(6) tran¬ 
sition. Since the E(5) symmetry corresponds to the second order shape phase transition 
with no shape coexistence, then no such competition exists as in the X(5) symmetry. 
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Figure 1: Comparison of energy level statistics of the states = O’*" in E(5) symmetry 
with different numbers of energy levels: (a) for = 20, = 20, (b) for = 40, 

= 40^ and (c) for = 80, = 80. In all figures, the solid lines and dashed lines 

describe the GOE and Poisson statistics, respectively. 
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Figure 2: Comparison of energy level statistics of the states J'^ = 8^ and iF = 6 in X(5) 
symmetry with the same manually truncated quantum number = 30 and = 30, 
but different relative strength A/B\ (a) for A/B = 20, (b) for A/B = 30, (c) for A/B = 
60, and (d) for A/B = 100. 
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